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Abstract. In this paper, we construct Hamilton-Jacobi equations for a great variety of mechanical 
systems (nonholonomic systems subjected to linear or affinc constraints, dissipative systems sub- 
jected to external forces, time-dependent mechanical systems...). We recover all these, in principle, 
different cases using a unified framework based on skew-symmetric algebroids with a distinguished 
1-cocycle. Several examples illustrate the theory. 
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1. Introduction 

A fundamental requirement for new developments in mechanics is to unravel the geometry that 
underlies different dynamical systems, especially mechanical systems. There are several reasons why 
this geometrical understanding is fundamental. First, it is a key tool for reduction by symmetries and 



2000 Mathematics Subject Classification. 70H20,70G45,70F25,37J60. 

Key words and phrases, skew-symmetric algebroids, 1-cocycle, linear Poisson structures, AV-bundlcs, Hamiltonian 
dynamics, Hamilton-Jacobi equation, affine nonholonomic systems, linear external forces. 

This work has been partially supported by MEC(Spain) Grants MTM2006-03322, MTM2007-62478, MTM2009- 
13383, ACIISI SOLSUBC200801000238, project IngenioMathematica(i-MATH) No. CSD2006-00032 (Consolider- 
Ingenio2010) and S-0505/ESP/0158 of the CAM. P.B. thanks CNPq(Brazil) for financial support and the Centre 
Interfacultaire Bernoulli (Switzerland) for its hospitality during the Program Advances in the Theory of Control Signals 
and Systems with Physical Modeling, where part of this work was done. The authors wish to thank Jose Mendez for 
helpful comments about the examples contained in this paper. 

1 



2 



P. BALSEIRO, J. C. MARRERO, D. MARTlN DE DIEGO, AND E. PADRON 



for the geometric characterization of the integrability and stability theories. Second, the effective use 
of numerical techniques is often based on the comprehension of the fundamental structures appearing 
in the dynamics of mechanical and control systems. In fact, the geometrical analysis of such systems 
reveals what they have in common and indicates the most suitable strategy to analyze their solutions. 
Finally, the geometrical approach has provided substantial contributions to neighboring areas, such 
as molecular systems, classical field theories, control theory, engineering, etc. 

Recent efforts have led to a unified framework for geometric mechanics based on a new structure, 
namely a Lie algebroid (see Section [2j), which represents the phase space for lagrangian mechanics 
and whose dual is the phase space for hamiltonian mechanics. These ideas were introduced in a 
pioneering paper by A. Weinstein |41) (see also [25]) where the equations of motion were derived from 
a Lagrangian function given on a Lie algebroid. This was done using the linear Poisson structure on 
the dual of the Lie algebroid and the Legendre transformation associated with that regular Lagrangian. 
The unifying feature of the Lie-algebroid formalism is particularly relevant for the class of Lagrangian 
systems invariant under the action of a Lie group of symmetries (see 23 for a survey on the subject; 
see also [3 ESI). 

As it turns out, the Lie-algebroid scheme is not general enough to include some interesting me- 
chanical systems. On a Lie algebroid, the Jacobi identity for the bracket of sections implies the 
preservation of the associated linear Poisson bracket on its dual. However, many interesting examples 
are not covered by this strong assumption, for instance nonholonomic mechanics (see HU EQl HQ] 
and references therein). Moreover, it would be interesting to find a general setting encompassing also 
some cases of dissipation of energy (for instance, explicit time-dependent systems, systems subjected 
to external forces or mechanical systems subjected to affine nonholonomic constraints). These reasons 
are our main motivation for introducing hamiltonian mechanics on more general objects, namely skew- 
symmetric algebroids equipped with a 1-cocycle; skew-symmetric algebroids will allow us to avoid the 
preservation of the Poisson bracket [3J [TOl HH [T3J, EH [37J and the 1-cocycle will introduce a dissipative 
character to the dynamics (for the geometric description of time-dependent mechanics, in terms of 
Lie affgebroids or, equivalently, in terms of Lie algebroids with a 1-cocycle, see [TT1 [T51 [T^l |2"71 12"51 1317] ; 
see also [3S]). Other approaches to the study of nonholonomic mechanical systems subjected to linear 
constraints, in the algebroid setting, have been also discussed in some recent papers (see [8l 19] 1311 132] ) . 
In these papers, the key tool is the notion of the prolongation of a Lie algebroid over a fibration. 

Our main goal is to derive a Hamilton-Jacobi equation for the case of skew-symmetric algebroids 
with a 1-cocycle. As it is well known, Hamilton-Jacobi theory for unconstrained systems is a useful 
tool for the exact integration of Hamilton's equations, for instance using the technique of separation 
of variables (see [I] and references therein). In other cases, this theory allows us to simplify the 
integration of Hamilton's equations or, at least, to find some particular solutions of the system. To 
summarize the idea for classical Hamilton's equations, consider a configuration manifold Q and a 
hamiltonian H : T*Q — > K. Then the Hamilton-Jacobi equation can be written as 

■"(?) t; — ) = constant 
oq 

for some function W : Q — > R. If we find such a function W, then the integration of Hamilton's 
equations (for initial condition along dW(Q)) is reduced to knowing the integral curves of a vector 
field on Q, defined as X H W = Tt t *q o X h o dW e 3£(Q), where t t »q : T*Q — >• Q is the canonical 
projection and Xh is the hamiltonian vector field associated to H . Hence, from the integration of a 
vector field on the configuration space it is possible to recover some of the solutions of the original 
hamiltonian system. 

A similar idea is also present in riemannian geometry when we look for a vector field X G 3L(Q) 
verifying V^A = (a geodesic or auto-parallel vector field), where V s is the Levi-Civita connection 
associated to a riemannian metric S on Q. Their integral curves are geodesies, that is, solutions of the 
geodesic second-order differential equations corresponding to S with initial conditions on Im X C TQ. 
Observe that, in general, X is not the gradient with respect to S of a function W € C°°(Q), which 
would be the case if we applied the classical Hamilton-Jacobi theorem. Hence, to recover this situation 
it is necessary to generalize the classical Hamilton-Jacobi equations. 
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On the other hand, recently, some of the authors of this paper proposed a generalization of 
the Hamilton- Jacobi equation for skew-symmetric algebroids (see [22]). Roughly speaking, a skew- 
symmetric algebroid is a vector bundle te ■ E — > Q equipped with a skew-symmetric bilinear bracket 
of sections and a vector bundle morphism, pe : E — s- TQ (the anchor map), satisfying a Leibniz-type 
property, i.e., a Lie algebroid structure without the integrability property, (for more details, see Sec- 
tion [2| . The existence of such a structure on E is equivalent to the existence of a linear almost 
Poisson bracket on the dual bundle te* '■ E* — > Q, or the existence of an almost differential d E on 
te ■ D — > Q which satisfies all the properties of an standard differential except that (d E ) 2 is not, in 
general, zero. 

Skew-symmetric algebroids were used in |22| to describe the Hamilton- Jacobi equation of nonholo- 
nomic mechanical systems. In this case E is determined by the linear constraints and the function 
W £ C°°(Q) is replaced by a 1-cocycle on the dual bundle E* (i.e., a section a of E* such that 
d E a — 0). With these ideas, one derives a Hamilton- Jacobi equation for nonholonomic dynamics, 
illustrating the utility of this new theory for the integration of different nonholonomic problems. 
Hamilton-Jacobi theory for standard nonholonomic mechanical systems has been also discussed in 
recent papers (see [6] HT] 136]). 

In this paper, we develop a Hamilton-Jacobi theory including, as particular cases, the Hamilton- 
Jacobi equation for skew-symmetric algebroids introduced in |22| and the case of auto-parallel vec- 
tor fields in riemannian geometry (Example 4.3 1, as well as a great variety of new examples (time- 
dependent hamiltonian systems, systems with external forces, nonholonomic mechanics with affine 
constraints...). With this objective in mind, we obtain the main result of our paper, Theorem 3.1 the 
Hamilton-Jacobi equation for a hamiltonian system on a skew-symmetric algebroid with a 1-cocycle. 
Moreover, our construction is preserved under the natural morphisms of the theory. This fact is 
proved in Theorem |3.6| We remark that this new version of the Hamilton-Jacobi equation is much 
more general than the one developed in |22| . since here, we do not require the 1-section solutions 
of the Hamilton-Jacobi equation to be closed. This fact is extensively used in Example |4.12| where 
we find solutions for the problem of a rolling ball in a rotating plane with time-dependent angular 
velocity looking for functions W £ C°°(Q) which do not satisfy (d E ) 2 W — (and thus, out of the 
cases studied in [22]). Moreover, the proof of the Theorem 3.1 is simpler and completely independent 
of the one done in [22] . 

The paper is structured as follows. In Section [2] we discuss some aspects of the geometry of 
skew-symmetric algebroids in the presence of a 1-cocycle. Moreover, given a hamiltonian section h of 
the AV-bundle associated with the skew-symmetric algebroid and the 1-cocycle, we obtain Hamilton 
equations for h. In Section[3] we formulate and prove the Hamilton-Jacobi Theorem for a hamiltonian 
system on a skew-symmetric algebroid with a 1-cocycle. In addition, we see that the Hamilton-Jacobi 
equation is preserved under the natural morphisms of the theory. Finally, in the last section, some 
theoretical and practical examples will illustrate the power of these new techniques as for instance: 
Hamilton-Jacobi equation for a particle on a vertical cylinder in a uniform gravitational field with 
friction, for a homogeneous rolling ball without sliding on a rotating table with time-dependent angular 
velocity or for the vertical rolling disk with external forces. 



2. Skew-symmetric algebroids, 1-cocycles and Hamiltonian dynamics 

2.1. Skew-symmetric algebroids and 1-cocycles. Let te ■ E — > Q be a vector bundle of rank 
n over the manifold Q. We denote by T(E) the C°°(Q)-module of sections of E. A skew- symmetric 
algebroid structure on E is a pair ([•, -],p), where [•, •] : T(E) x T(E) — > T(E) is a M-bilinear skew- 
symmetric bracket on T(E) and p : E — > TQ is a vector bundle map (the anchor map) such that 

h/7l = /h7]+pW(/)7, <t,^£T(E), f£C°°(Q). 
Note that p : E —> TQ induces a homomorphism of C 00 (Q)-modules which we denote also by p : 

T{E) -> X(Q) (see [3 m El ELU 133). 
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If the bracket [•, •] satisfies the Jacobi identity then (E, [•, •], p) is a Lie algebroid (see, for instance, 
26 ). In such a case, we have that the anchor map is a morphism of Lie algebras, i.e. 

p(l<r,i\) = \p(<r),p(rf)], for<7 )7 er(£). 

On a skew-symmetric algebroid structure ([•,•],/£>) on the vector bundle te ■ E — > Q it is induced 
the almost differential d E : T(A'E*) -> r(A ,+1 £'*) as a M-linear map given by 

(d E f)(a) = p(a)(f), 

(2.1) 

(d E a)(a ll ) = p(o-)(a(7))-p(7)(a(«r))-a([a-,7]) 

and 

d E (fa A fa) = d E fa A fa + {-if fa A d E fa, 

for / e c*°°(g), a e r(£*), cr, 7 e t(e), fa e r{A k E*) and fa e t(a'E*). 

Note that d E is defined in a similar way that the standard differential over a manifold. However, 
there are important differences between them. In what follows, we will discuss some facts related with 
these differences. 

Firstly, unlike the case of the stardard differential on a manifold, we have that the almost differential 
d E of a skew-symmetric algebroid (E, [•, •], p) is not, in general, a cohomology operator, i.e., (d E ) 2 ^ 0. 
In fact, (d E ) 2 — if and only if ([•, •],/>) is a Lie algebroid structure. 

For the particular case of a function g £ C°°(Q), we deduce from (2.1) that (d E ) 2 g = if and only 
if X(g) — for all X £ D, where D is the finitely generated distribution given by 

D:= S pan{[p(a),p(j)}-p(la, 7 j) : a, 7 e T(E)} C X(Q). 

On the other hand, if Q is a connected manifold, in general, d E g = does not imply that g : Q —> M 
is constant. In other words, if Q is a connected manifold, in general, the vector space 

H°(d E ) = {/ G C*°°(Q) such that d £ / = 0}, 

is not isomorphic to K. Note that, when E is a Lie algebroid, H°(d E ) is the Lie algebroid cohomology 
0-group of E. Even in this case, one can not guarantee that H°(d E ) is isomorphic K. However, if Q 
is connected and J5 is transitive, i.e., p{E) — TQ, then H°(d E ) = K. 

In [22] the authors discuss the relation between a function <? £ C°°(Q) being constant and d E g = 0. 
In order to remember these results we introduce the notion of completely nonholonomic distribution 
(see [33]). 

Definition 2.1. A distribution D C TQ is called completely nonholonomic (or bracket generating) 
if {Xk, [Xk,Xi], [Xi, [Xk,Xi]], ... £ X(Q) : Xj(q) £ D q \/q £ Q} spans the tangent bundle TQ. 

The Lie brackets of vectors fields in D generate a flag D C D 2 C ... C TQ with 

T> 2 = D + [T>,T>], T) r+1 = T> r + [D, T> r ] 

where 

[D, D fe ] = span{[X,F] : X £ D and Y £ D k } 

with the spans taken over smooth functions on Q (for details, sec 33J ) . 

Here, we have two extreme cases: on one hand, the distribution D can be involutive, then we have 
D = D 2 = D r , W £ N>o- On the other hand, if D is completely nonholonomic, then there exists 
r £ N >0 such that D r = TQ. 

Now, consider a skew-symmetric algebroid (E, [•, -],p) on a manifold Q and the following finitely 
generated distribution T> given by 

D q =p q (E q ), for all q eQ. 

Proposition 2.2. [22j // (iJ, [•,•], p) is a skew- symmetric algebroid over a connected manifold Q, 
such that D = p(E) C TQ is a completely nonholonomic distribution, then H°(d E ) = K. 
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However, in some examples, the distribution D is not completely nonholonomic. In such a case 
there is r 6 N >0 such that U r - 1 C D r = U r+1 C TQ. This distribution D r is the smallest Lie 
subalgebra of X(Q) containing 2). Let us consider the associated generalized foliation over Q. The 
leaf of this foliation over a point q E Q, is just the orbit 

L = {(j)?* o ...ocj)fi(q) eQ : f,€l, andl.eD with i = l,...,k, k G N >0 } 
where </>^ is the flow of the vector field Xi at time ti(see [2"].|39]). 

Theorem 2.3. |22| Let (E, [•, 6e a skew- symmetric algebroid over a manifold Q. Consider the 
leaf L of D r as described above. Then 

(i) It is induced a skew- symmetric algebroid structure ([■, Pl) on the vector bundle Tl '■ El — > 
L with E^ := U qe LE q . Moreover, the distribution Pl{El) on L is completely nonholonomic. 

(ii) If f E C°°(Q) is a function such that d E f = 0, then its restriction to L is constant. 

Next, we will see that any skew-symmetric algebroid structure ([•, -J, p) on the vector bundle te '■ 
E — > Q induces an almost Poisson linear bracket {•, •} : C°°(E*) x C°°(E*) — > C°°(E*) on the space 
of functions on E* , that is, {•, •} is a skew-symmetric M-bilinear bracket which is a derivation in each 
argument with respect to the standard product of functions and with the extra property that the 
bracket of two linear functions is again a linear function. Indeed, this bracket is characterized by the 
following relations (see [3 Q2 QH [HI [22] ) 



{ct,7} = -[cr,7l, {foTE^a} = p(a)(f) o t e . 

(2.2) 

{f °T E ,,goT E *} = 0, 

for all cr, 7 € T(E) and f,g € C°°(Q) and where £ : E* —> M is the linear function associated with the 
section ( eT(E). 

Now, we will endow our skew-symmetric algebroid with an additional structure: a distinguished 
section (p £ T(E*) which allows us to consider some interesting examples. 

Let us consider a section <f> of E* . Denote by <fi v € X(E*) the vertical lift of the section <j> £ T(E*), 
that is, the vector field defined by 

^ (a) = (<t>(T B ,(a))) v a , VaeE*, 

where ^ : E*^,^ — > T a (E* ^ ) is the canonical isomorphism between the vector spaces E* B ,r a ) 
and T a (E;^ (q) ) . Note that 

v (ct) = <j)(a) o t e * = a o 4>o t e *, (2.3) 

for all a G r(S). 



Using (2.1), (2.2) and ( |2.3[ ), we obtain the following formula which describes the differential of 
in terms of the linear bracket {•, •} on E* 

d s ^( 7 ,a)or £ , = -({0 v ( 7 ),a} + { 7 ,0 v (a)}-0 v ({7,a}) 

(2.4) 

= -({7) O- o o T E * } - {7, d} + {7 o o r B « , cr}), 

for all 7, cr G r(E). Thus, 4> £ r(S*) is a 1-cocycle, i.e., d E 4> = 0, if and only if 

^({^1^2}) = {4> V (<Pi), ^2} + W l7 ^( V2 )}, (2.5) 
for all (fi,f2 G C°°(E*). Moreover, equation (|2.5[) is equivalent to the fact that the linear bivector 



n^* on E* associated with the bracket {•, •} is invariant with respect to the vector field 4> y G X(E*), 



i.e., 

V n £*=0- (2-6) 
In fact, from (2.2), (2.5) and (2.6) one may conclude the following result: 

Proposition 2.4. Let te : E — > Q be a vector bundle over Q and </> G T(E*) a section of the dual 
bundle of E. Then, the following statements are equivalent: 

(i) E admits a skew-symmetric algebroid structure ([•, -],p) such that d E <fi = 0. 
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(ii) There is a linear bivector He* on E* which is invariant with respect to the vertical lift 

f el(£*)»/f 

2.2. Hamiltonian dynamics. Let (E, [•, •],/?) be a skew-symmetric algebroid over Q of rank n, and 
cf> £ T(E*) be a 1-cocycle of E. Denote by <fr : E — > E the corresponding linear function induced by <fi 
on E and suppose that, for all q £ Q, <p\E q 7^ 0. Then, one may consider the affine bundle 

ta:A:= <£ _1 (1) -> Q 

of rank n — 1 with associated vector bundle Ty : V := <j)~ 1 (0) — > Q. Note that V is a skew-symmetric 
algebroid over Q with structure ([•, -]y, py) given by 

«v o [ct,7]v = [«y(c r ),^(7)] I /fy(o-) = p(*v(o")) 

for all cr, 7 G r(V), where iy : V" —> E is the canonical inclusion. Thus, we have the corresponding 
linear almost Poisson 2- vector IIy» on V* . 

On the other hand, the map \i := i* v : E* — > V* defines an affine bundle of rank 1 modeled over 
the trivial vector bundle pr\ : V* x E — >• V* (an AV-bundle in the terminology of [H]). Using (2.2 1 
and the fact that the canonical inclusion iy : V — > E is a skew-symmetric algebroid monomorphism, 
we deduce that \i : E* — > V* is an almost Poisson morphism. Thus, if a q £ E* ql we have that the 
following diagram is commutative 



T* E* 



T, 



T a E* 
V* 



Here, #n E , : T*E* ^ TE* (respectively, : T*l/* TU*) is the morphism of C°°(Q)- 

modules induced by the almost Poisson bivector n^;. (respectively, Hy*). 

Using again (2.2), we also deduce that 

d v f = fiod E f &ndd v (noa) = A 2 fiod E a, f £ C°°(M), a £ T(E*), 

where A 2 /i : A 2 E* — > A 2 V* is the extension of /z to the corresponding vector bundles. 

Note that the set of the global sections r(/z), of fx, is an affine space modeled over C°°(V*). In 
addition, if h £ then fJ,(a q — h(fx(a q ))) — 0, for q £ Q and a q £ E*. Thus, one may define a 

function Fh £ C°°(E*) characterized by 



a q - h{j±(a q )) = F h (a q )cj>{q), 



(2.7) 



for all q £ Q and for all a q £ E* Moreover, we have 



(4> V (F h Ma q ) 



dt\t=o 



F h (a q + tq>(q)) = 1. 



Therefore, it follows that (j> v (F^) — 1. In fact, there is a one-to-one correspondence between r(/x) and 
the set of functions F on E* which satisfy the relation 



(2.8) 



(see HU). 



In what follows, we will associate to each section h £ r(/z), a vector field on V* . From (2.5) and 



(2.8), we deduce that, for every section h of the bundle fx : E* — > V* and G £ C°°(V*), the function 
{Go n,Fh} is /i-projectable (note that KerT aq n —< <f> v (a q ) >). Thus, for each h £ T(/i) we can 
consider a vector field Rh on V* which is characterized by 



R h (G)o f , = {Go f x,F h }, G£C°°(V*)- 



(2.9) 
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This vector field is called the hamiltonian vector field associated with the section h. If !K^ E * E X(E*) 
is the hamiltonian vector field associated with the function Fh <E C°°(E*) with respect to the linear 
almost Poisson bracket He* , i.e., 

then, from (2.9), we deduce that 

R h op = TpoOi^'. (2.10) 
The integral curves of the vector field Rh are the solutions of the Hamilton equations for h. 

2.3. Local expressions. Let E be a vector bundle on Q of rank n, with a skew-symmetric algebroid 
structure ([•,•], p). 

Fixed a section <fi of E* such that d E (f> — and 4>\e f° r a U <7 € Q. Then it is induced a local 
basis {eo, ea}a=i,...,n-i of E adapted to the 1- section in the sense that <f>(eo) — 1 and <p(e a ) = 0. In 
terms of this basis we have the local structure functions, Q c ab , p l a , Cqj,, p l € C°°(Q) of E defined by 

9 9 
\e a ,e b \ = e c ab e c , [e ,e b ] = e c ob e c and p(e a ) = Pag^: P( e o) = Po^- 

Note that the condition cP0 = implies that C" fc = 0, for all a, b £ {0, . . . , n — 1}. 

Moreover, with respect to the induced local coordinates (q 1 ,po,p a ) on E* , the local expressions of 
the vector field </> v € X(E*) and the linear almost Poisson bivector He* are 

9po 

T-r i 9 3 , 9 9 „ 9 9 1 „ 9 9 

"dq l dp dq l dp a u dp op b 2 dp a dp b 

If (q l ,p a ) are the corresponding coordinates of V* , the local expression of p : E* — > V* is 

p{q\po,p a ) = {q\Pa)- 

Let h : V* —> E* be a section of whose local expression is 

h{q\p a ) = {q\-H{q\p b ),p a ) 
where H is a local function of V* . Then the associated function F b : E* — > K is 

F h {q\p ,p a ) = Po + H(q\ Pa ). (2.11) 

Moreover, the local expression of the Hamiltonian vector field associated with this section h : V* — > 
E* is given by 

R h - (po + Pl w ) w + (-Pl w + (e& + ei b — ) Pc )-. 

Thus, the Hamilton equations are 

dq l , , dH dp b ;dH nc dH . 

7ft = * + ^ = "'W + (CSb + ^dp-J^ 

A Lagrangian version of these equations was considered in [TB] (see also [55]). 

It is important to note that the previous dynamics on V* has a dissipative character. In fact, in 
the case when the AV-bundlc p : E* — > V* is trivial, then the local function H is global and it is the 
hamiltonian function on V*. In addition, 

R h {H) o p = {H o p, F h } ^ 0. 

The local expression of this dissipative term is 

OH OH 
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2.4. Examples. Next, we will describe two interesting examples of skew-symmetric algebroids which 
will be useful for the mathematical description of the mechanical systems considered in this paper. 

Example 2.5. Consider a Lie algebroid structure ([•,•], p) (or more generally a skew-symmetric 
algebroid structure) on a vector bundle Tg : E — > Q of rank n — 1 and F : E — > E a homomorphism of 
vector bundles (over the identity of Q). Then, on the vector bundle t Rx ^ : R x E —> Q, it is induced 
a skew-symmetric algebroid structure (E :=R x E, [•, -IrxEtPrxe) given by 

[(/,*), (0, 7)]»xis = (p(v)(9)-p(7)(fU*,7l+9F(<r)-fFm, 

PrxeH^) = Pi' 7 ), 

for all (f,tr), (g,j) e T(E) = C°°(Q) x T(E) = T(R x E). 

Note that the space T(A 2 (IR X_E*)) may be_ identified with T(E*) © T(A 2 E*). Under this identifi- 
cation, for (/, a) £ C°°(Q) x r(i?*) = r(Ix£'), we obtain that 

cF xB (f ia ) = (F*a~d s f,d E a), (2.13) 

where F*a is a section of £7* defined by 

(F*a)(a) = a(F(a)), Vct € T(E). 



From (2.13), we have that (1,0) £ r(R x E*) is a 1-cocycle, and its vertical lift is just the vector 
field on R x with po the global coordinate on R. 

In this case, the linear almost Poisson bivector on R x E* is given by 

n ix£« = n E* + A ( F *) V > 

where IT^, is the linear Poisson bivector on E* induced by the Lie algebroid structure on E and (F*) v 
is the vector field on E* defined by 

(F*) v : E* -> TE* , (F*Y (a) = F*(a) v a . 
Here, F* : E* — >• is the dual morphism of F : E -> £7. 

On the other hand, the affine bundle associated with (R x E, [•, ■] RxJ j, /% x .e> (1j 0)) is just ^ : 
{1} x E = E ^ Q with associated vector bundle T{ } x ^ : {0} x E = £ — » Q. 

Moreover, the associated AV-bundle p : R x £7* — > E* is just the trivial affine bundle over E* of 
rank 1. Therefore, there is a one-to-one correspondence between sections h : E* — > R x E* of this 
affine bundle and functions H : E* — > R. Furthermore, the hamiltonian vector field associated with ft, 
is just 

R h = 'K^' -(F*Y, (2.14) 

where J{^ E * G 3t(E*) is the hamiltonian vector field associated with H with respect to the linear 
Poisson structure il^, . 

Note that, in this case, the dissipative term is given by 

R h {H) = -(F*Y (H). 

Now, we will give some expressions in coordinates of the above objects. Let us consider local 
coordinates (q l ) on the manifold Q and a local basis {e a } a =i,...,n-i of T(E). 

A local basis of sections of E — R x E is {eo, e a } a =i,...,n-i, where eo = (1,0) and e a = (0, e Q ). 
Locally, the homomorphism of vector bundles F : E — > E is given by the functions F^ £ C°°(Q), 
where F(e a ) = F%e b . Then, 

[eo,e a ] RX £; = — F a ei,, [ e a,e&] RX £ = C^ee 

PE{e ) = 0, Pij(ea) = p\-§-[. 

The linear almost Poisson bivector on E* = R x £7* is now: 

„ ,9 9 9 9 1 d d 

9<f dp a dpo 9p 6 2 dp a 9p 6 
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f) 



Given a hamiltonian function H : E* — > K then the Hamilton equations of motion are: 

dcf_ i dH_ 

dp b i^H dH 

dF = ~ pl W + GabPc dp- a ~ FbPa 

Example 2.6. (see [22]) Let (E, [•, •], p) be a Lie algebroid on a manifold Q. Suppose that we have a 
vector subbundle D of E and a projector, i.e., a vector bundle morphism O 3 : E — > D (over the identity 
of Q) such that 7\ D = id. Denote by io '■ D — > E the natural inclusion from D to E. Then, we may 
induce a skew-symmetric algebroid structure on D as follows 

{v,1\d = ^(Pd oa,i D o 7]), p D (a) = p(i D o a), 

for all <7, 7 € r(D). Note that, in general, ([•, -]d,Pd) is not a Lie algebroid structure on D. 



3. Hamilton- Jacobi equation, skew-symmetric algebroids with a 1-cocycle and linear 

almost poisson morphisms 

3.1. Hamilton- Jacobi equation. Let te : E —> Q be a vector bundle, of rank n, on the manifold 
Q with a skew-symmetric algebroid structure ([•, 

Consider e r(£'*) a 1-cocycle of -E 1 which satisfies the following condition: 

(j)\E q ^ 0, for all qeQ. 

Then, as we have shown in Section 



2.2 



the vector bundle ry : V := <j) 1 (0) — > Q of rank n — 1 admits 
a skew-symmetric algebroid structure which we denote by ([•, -]y, pv)- 

If h is a section of the corresponding AV-bundlc p : E* — > V* then (E, [•, •], p, <j>, h) is said to be a 
Hamiltonian system on E . 

In such a case, for each section a of V* , one may define a vector field on Q as follows 

K = Try, oR h oa, (3.1) 

where Rh is the hamiltonian vector field associated with the section h. 

On the other hand, we may introduce the following map 

~ : ^(E*) -> T(E), w^w, 

where ui is characterized by 

/3{u) =uj(J3 v )ohoa, for all /3 e T(E*). (3.2) 

We remark that 2 € r(£) since, if / € C°°(Q) then (//3) v = (/ o t e *)/3 v and 

T£;» o h o a = id. (3-3) 

Moreover, it follows that 

dff = 7, d(foTE*) = 0, and Fw = (f o /1 o (3.4) 

for 7 G r(^), / G C°°(Q), F G C°°(E*) and w G O 1 ^*)- We wiU denote by ^ the section of E 
given by 

C - dF h . (3.5) 

Now, we state the main result of this paper. 

Theorem 3.1. Let (E. [•,•], p, <fi, h) be a Hamiltonian system on E. If a £ T(V*), we have that the 
following statements are equivalent: 

(i) If c : I Q is an integral curve of i?^ G X(Q) then a o c : I — >• V* is an integral curve of 
Rh G X{V*). 

(ii) a G T(V*) satisfies the Hamilton-Jacobi equation, i.e., 

po(i C c.d E {hoa)) =0. 
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Proof. From fl2.2| ), ( |2.4[ ) and (|3_3|), we deduce that for all 7, a e T(E) 

d E (h o a) (7, a) = (p(j)(& h o a) o te* + {7, <?} — p{cr)^y ° h o a) o te*) ° ft ° a 

=— n^* (cfy, (ho a o te*)* (da)) o h o a + Ue* (d/y - (ft o a o t^* )*(<fy), d£r) o ho a. 

Moreover, using ( |3.4[ ), it follows that 

= d E (hoa)(d(fo^7 E *)idF) 



and by (2.2 1, (3.3) we obtain also that 

= -IlE*(d(foTE*-)AhoaoTE*)*(dF))ohoa 

+U E *(d(f ot e *) - (hoao T E *)*(d(f oTE*)),dF) oho a 
for all F G C°°(E*) and / G C°°(Q). 



(3.6) 

(3.7) 
(3.8) 



Therefore, from (3.3), (3.4), (3.6), (3.7) and (3.8), we conclude that for all wi,w 2 G ft 1 ^*) 



d E (ho a)(uj 1 ,uj 2 ) 



-He* (wi, (ft. o a o te»)*W2) o ho a 
-He* (wi — (ft o a o te* )*wi, W2) o ho a. 



(3.9) 



Next, we will denote by iy : V — )• E the natural inclusion. Then, if cr G T(V), considering in (3.9) 
the particular case when o>i = dFh and U2 — d(p o /i) = e2(iy o cr), we have that 



d E (h o a)(Ch , iy o <t) = M F ^* (cr o a o ry» ° p) ° ho a 
Note that (ft o a o T E *)*(dF h ) = (see dO). 



(cr o (tx) o ft o a. 



Now, using (2.10) we conclude that 

d E (h o a)(Ch, i V oa) = (Ta o Rf L )(d) - R h (a) o a. 

Next, we remark that statement (i) in the theorem is equivalent to 

TaoRf^ R h o a. 



So, if this relation holds then, from (3.10), we deduce that /1 o (i^d E (h o a)) = 0. 
Conversely, suppose that 

d E (h o a)(Ch,iv o d) = 0, for all cr G T(V). 



(3.10) 
(3.11) 

(3.12) 



In order to prove (3.11 ), it is sufficient to see that the following two relations are satisfied: 

(a) R%(a oa) = R h (a) o a, for all cr G T(V), 

(b) R%(f) = R h (fo Tv *)oa, for all / G C°°(Q). 



Statement (a) is a consequence of Eq. (3.1) and (b) follows from (3.10) and (3.12). 



□ 



In what follows, we write the local expression of the Hamilton- Jacobi equation. 

Consider local coordinates (q l ) on Q and a local basis {eo,e a } a =i....,n-i of T(E) adapted to the 
1-cocycle <j> as in Subsection |2.3| 

Denote by (q l ,po,p a ) (respectively, (c/ l ,p a )) the corresponding local coordinates on E* (respectively, 
V*). Then, the section a : Q — > V* and the hamiltonian section ft : V* — > E* are written in terms of 
these coordinates as 

a(q l ) = (q\a a (q i )), h(q\p a ) = (q\ -H(q i ,p a ),p a ). 

On the other hand, a 1-form u G il 1 (i?*) can be writt en i n these coordinates as u = u) % dq % + 
uj a dp a + u°dpo with w\w a , w° G C°°(£'*). Therefore, from (Q we obtain that the section u> G r(J3) 
is given by 

cD = (uj a o ft o a) e a + (cj° o ft o a) eo- 

Thus, 



eo + ( q — ° a)e a . 

OPa 



(3.13) 
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Now, from (2.1) and (3.13), the Hamilton Jacobi equation is given locally as follows 

dH 



i i dH da a 



(Pa 



dq i 



OH 



for all a = 1, . . . , n — 1. 



To finish this subsection, we will show some consequences of Theorem |3.1| which will be useful for 
the next examples. 

Corollary 3.2. Let (E, [•, •], p, <f>, h) be a Hamiltonian system on E. If ft £ T(E*), then the following 
statements are equivalent: 

(i) If c : I — > Q is an integral curve of Rf°^ £ X(Q) then p o f3 o c : I V* is an integral curve 
of lh, ■ X:\-U 

(ii) (3 £ T(E*) satisfies the Hamilton- Jacobi equation, i.e., 

/i o i „of< d E f3 + d v (F h o(3)=Q. 



Proof. Using (2.7), we deduce that the Hamilton- Jacobi equation for fi o fj g r(V*) is 

pU/(MF k o^)) = 0. 



Since is a 1-cocycle then 



d E ((3 - (F h o p)c£) = d E [3 - d E (F h o f3) A <fi. 



(3.14) 
(3.15) 



On the other hand, from (3.2), it follows that 



i 



and therefore, using (3.15), we obtain that (3.14) is equivalent to 

/i o i^of,d E f3 + /io d E (F h ofi)- (i^op (d E {F h o p)))n o = 0. 
Finally, the corollary is an immediate consequence of Theorem |3.1| and the relations 



Hod E (F h 



d v {F h o(3) and fj, o = 0. 



□ 



Corollary 3.3. Let (E, [•, •], p, <f>, h) be a Hamiltonian system on the vector bundle te '■ E — > Q on 
the connected manifold Q. Suppose that the finitely generated distribution V defined by V q :— pv(Vq) 
for all q G Q, is a completely nonholonomic distribution. If /3 €E T(E*) is a 1-cocycle of E* , then the 
following statements are equivalent: 

(i) If c : I — > Q is an integral curve of G X (Q) then fi o f3 o c : I — > V* is an integral curve 

ofR h eX{V*). 
(ii) (3 £ T(E*) satisfies the Hamilton- Jacobi equation 

Fh ° — constant. 

3.2. Linear almost Poisson morphisms and Hamilton-Jacobi equation. As we pointed out 
in the introduction, one important advantage of dealing with unconstrained hamiltonian systems on 
Lie algebroids, or constrained systems on skew-symmetric algebroids, is that reduction by symmetries 
can be naturally handled by considering morphisms between Lie algebroids [H] (see also [7]), or 
respectively morphism between skew-symmetric algebroids, [22] . In the following section we deal with 
morphisms between skew-symmetric algebroids with a 1-cocycle, in order to show that the Hamilton 
Jacobi equation is preserved by such morphisms. 

Suppose that te ■ E — > Q and : E — > Q are vector bundles over Q and Q, respectively. Consider 
a vector bundle morphism (vt','0) between E* and E* 
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Denote by A $ : A k E* — > A k E* the corresponding vector bundle morphism on ip : Q — > Q, induced 
by the pair (^,ip), between the vector bundles A k E* — > Q and A k E* — > Q. A section a G T(A k E*) is 
($>,ip) -related with a section a € T(A k E*) if 

A fc \I> o a = a o ip. 

Definition 3.4. Let (E, [■, -Je, p, <p, h) and (E, [•, J E , p, <p, h) be hamiltonian systems on E and E, 
respectively. Suppose that (^,ip) is a vector bundle morphism between E* and E* . Then, the pair 
(&,ip) * s sa id to be a hamiltonian morphism if: 

(i) (\P, ip) is an almost Poisson morphism, that is, 

{Fio*,F 2 o*}p ={F U F 2 } E « o% for F 1 ,F 2 eC°°(E*) 

where {•, -}e* (respectively, {•, -}e*) * s the linear almost Poisson bracket on E* (respectively, 
E*); 

(ii) <fi and <p> are ip) — related and 

(hi) F- h om = F h . 

Now, we prove the following result 

Proposition 3.5. Let (E, [•, p, (p, h) and (E, [•, -J E , p, (p, h) be hamiltonian systems on E and E, 
respectively, and ($!,ip) be a hamiltonian morphism between E* and E* . Then: 



(i) There exists a linear almost Poisson morphism ^ : V* 
diagram is commutative 



E* 



E* 



>- v* 



V* ( over ip ) such that the following 

(3.16) 











L n 






Y 



v* 



(ii) The vector fields Rh G X(V*) and Rj. L G X(V*) are V?— related, that is, 

R- h o^ ^T^ o R h . 

(iii) If a € r(V*) and a G T(V*) are (\& , ip) -related then the vector fields R^ G X(Q) and 
Rj^ G X(Q) are ip-related, that is, 

Rf i0 yj = TyjoR%. 

Proof, (i) Using that ip) is a vector bundle morphism and the fact that ^ o (p = <j> o tp, i t follo ws 
that there exists a vector bundle morphism ^ : V* — > V* (over ip) such that the diagram (3.16) is 
commutative. Moreover, since "J 7 ,/! and p. are linear almost Poisson morphisms, we deduce that \& 
also is a linear almost Poisson morphism. 



(ii) The condition F^oty = implies that 



o v]/ 



Ph. 



(note that \t is an almost Poisson morphism). Thus, from (i) and (2.12), we have that 



R- h o * = T^!oR h . 



(iii) Using (i), (3.1 ) and the fact that \toa = ao^,we conclude that the vector fields R% and R? 
are ip-rel&ted. 



From Proposition |3.5[ we deduce that following result 



□ 



Theorem 3.6. Let (E, [•, -Jb, p, cp, h) and (E,l-,-J E ,p,(p,h) be hamiltonian systems on E and E, 
respectively, and (^,ip) be a hamiltonian morphism between E* and E* . If a G T(V*) satisfies the 
Hamilton- J acobi equation for h, a G T(V*) is (\I/, ip) -related with a and ip is a surjective map then a 
satisfies the Hamilton- Jacobi equation for h. 
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Remark 3.7. Let (E, [•, -Je, p, cf>, h) and (E,l-,-J E ,p,4>,h) be hamiltonian systems on E and E, 
respectively, and ('J, ip) be a hamiltonian morphism between E* and E* . Suppose that a € T(E*) is a 
1-cocycle of E* such that F^o a — constant. Then, if a € T(E*) is a 1-cocycle of E* which is ip)- 
related with a, it is clear that F^oa — constant and, therefore, a is a solution of the Hamilton- Jacobi 
equation for h. o 

4. Examples 

In this section we will apply our theory to two type of mechanical systems: uncontrained mechanical 
systems with a dissipative character (with linear external forces or time-dependent systems) and 
nonholonomic mechanical systems subjected to affine constraints. In the last part of the section we 
will discuss the case of a nonholonomic mechanical system with external linear forces. 

4.1. Uncontrained mechanical systems with a dissipative character. 

4.1.1. Mechanical systems on Lie algebroids with linear external forces. (See Example |2.5[ ). Let us 
consider a Lie algebroid structure (or more generally a skew-symmetric algebroid) ([•,•], p) on a vector 
bundle t e : E ~r Q and a homomorphism of vector bundles F : E — > E (over the identity of Q). With 
these ingredients, it is induced on the vector bundle t^ xE : M x E — > Q, a skew-symmetric algebroid 
structure (E := Rx E, [•, -1 Rx e,Prxe) such that (M) € r(Rx E*) = C°°(Af) xT(E*) is a 1-cocycle. 

Let H : E* — > R be a differentiable function (Hamiltonian function) on E* . Denote by h : E* — > 
Kx£* the induced section of = pr 2 : 1 x £* -) £* by H, i.e., 

/i()9,) = (-H((3 q ), fi q ), for all q G Q and q e £*. 



The vector field on is just J£^ s * - (F*) v (sec (2.14)). Moreover, 



RZ = Tr E , o o a. 

On the other hand, for each a G r(i?*) one may define a section, of .E as follows 

/?(G)=/3 v (#)°a, 

for /3 <E r(25*). Then, under the identification T(R x E) = C°°(Q) x T(E), (% is just (1,Ch)- Thus > 
using Corollary |3.2| we conclude the following result 

Corollary 4.1. Let (E, [•,•], p) be a Lie algebroid (or more generally a skew- symmetric algebroid) 
with Hamiltonian function H : E* — » R, and F : E ^ E be a vector bundle homomorphism. If 
a G r(£J*), £/ie following statements are equivalent: 

(i) If c : I — > Q is an integral curve ofTr E , o o a e X(Q) i/ien aoc:/->£* is an integral 
curve o/5<£** - (F*) v . 

(m) a G r(25*) satisfies the Hamilton- Jacobi equation, i.e., 

i Cfi d S a + d S (Hoa)+ F*(a) = 0. 

Remark 4.2. (i) When a is a 1-cocycle and F = 0, we recover the result of [22 . Applications of 

this result to nonholonomic mechanical systems subjected to linear constraints were discussed 
there. Note that in this case the dissipative term is zero. 

(ii) In the particular case when E is the standard Lie algebroid ttq : TQ — > Q, then, using well- 
known results (see, for instance, |24j). we deduce that there exists a one-to-one correspondence 
between the vector bundle morphisms F : TQ — » T*Q (over the identity of Q) and the semi- 
basic 1-forms on TQ which are homogeneous of degree 1. A semi-basic 1-form j3 : TQ — > 
T*(TQ) on TQ is said to be homogeneous of degree 1 if £a/3 = P, where A is the Liouville 
vector field on TQ. Indeed, if 

F{q\q i ) = {q\F i j {q)cp) 
then the corresponding 1-form P on TQ is given by 
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Using this result and Corollary 4.1 we deduce Theorem 3.4 in [T7] for the particular case 



when the semi-basic 1-form f3 on TQ is homogeneous of degree 1. 

o 

Example 4.3. Standard mechanical systems. Let E be the standard Lie algebroid t^q : TQ — > 
Q. In this case the differential d E is the standard differential, d, on Q. Suppose that F = and that 
H : T*Q — > K is a hamiltonian function. If a is a 1-form on Q then the Hamilton- Jacobi equation is 

ixfjda + d(H o a) = 0, (4.1) 

where X% is the vector field on Q defined by X%{fi) = /3 V (H) o a, for all (3 G f2 1 (Q). 

If Q is connected and S : Q — > K is a function on Q, using the 1-form a = dS 1 , we obtain the 
standard Hamilton- Jacobi equation on Q, i.e., 

-ff o dS 1 = constant. 

On the other hand, let S be a riemannian metric on a n-dimensional manifold Q, i.e, a positive- 
definite symmetric (0, 2)-tensor on Q. The metric 9 induces the musical isomorphisms 

b s : X(Q) — > n\Q), \? 9 (X)(Y) = 9(X,Y), 
ttg :^(Q) — >3C(Q), « g (a) = bg 1 (a) 

where 1,^6 X(Q) and a G O^Q). 

Associated with the metric 9 there is an affine connection V s , the Levi-Civita connection, 
determined by: 

[X, Y] = V|Y - Vf-X (symmetry) 
X(9(Y, Z)) = 9(V| Y, Z) + 9(F, V|Z) (metricity) , 
for every X,Y,Z e X(Q). 

Considering a vector field X G X(Q) and the associated 1-form a — bg(X), we will analyze the 
meaning of the Hamilton-Jacobi Equation (4.1 ) for the Hamiltonian H : T*Q — > R defined by H(rj q ) 



^9*{r]q,T} q ), where 9* is the induced metric on T*Q and rj q G T*Q. First, we observe that the section 
Q of TQ is just the vector field X. Then, for Y G X{Q) 

da(X,Y)+Y(Hoa) - d(b s (X))(X, F) + ±Y(S(X,X)) 

- X(9(X, F)) - ±Y{S{X, X)) - 5(X, [X, Y]) 
= 9(V|X,F). 



Therefore, the Hamilton-Jacobi equation ( |4.1[ ) for the case of a Hamiltonian defined by a riemma- 
nian metric is equivalent to the condition for auto-parallelism of vector fields, that is, vector fields 
X G £(Q) such that V|-X = 0. 

Thus, if we have a vector field X which satisfies the auto-parallelism condition, each integral curve 
c : I — » Q (which is a geodesic) induces a solution of the Hamilton equations of the mechanical system, 
which is just 

\> 9 (X)oc: I ^T*Q. 

Example 4.4. The test particle under the gravitational interaction of two masses. Con- 
sider the problem of the motion of a particle moving under the gravitational effect of two masses 
ttl\ and 7722 , which in turn move in circular orbits about their common center of mass and are not 
influenced by the motion of the particle (classical planar circular restricted three-body problem). Take 
a coordinate system rotating about the common center of mass with the same frequency as the two 
masses so that both of them lie on the x-axis with coordinates (— /i2, 0) and (/j,i, 0), where fa — m ^ m2 
(see [ID OH]). 

The system is described by the Lagrangian function: 

L{x, y, x, y) = \{x - y) 2 + \(y + x) 2 - — - — , 
2 2 n r% 
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where r\ = (x + ^) 2 + y 2 and v\ = {x — /ii) 2 + y 2 . 

The equations of motion adding a drag force F = (F±(x, y, A, y), F 2 {x, y, A, y)) are (see |34j ) : 

au - 

x-2y-x = — — F\ , 

dU ~ 

y + 2x-y = -— F 2 , 

dy 

where U(x, y) = ^ + ^ and F : TR 2 -> T*R 2 . 

Now, we will describe this system in our geometric framework. Consider the homomorphism F : 
TR 2 -> TR 2 given by 



where F h a G C°°(IR 2 ). Then, on the vector bundle r : R x TR 2 -> R 2 it is induced a (transitive) skew 

I / - 'iliviinvr.u! ^'1 vn.'l lii'n /l/^, 'I'll wi/ ! ill i 1 i / i I / w • ■ I 1 I i - 1 ^' i ^ ,f / ■ ■ I I I M / ■ ■ ■ / I ) ___ j , . > , , j ^ 



symmetric algebroid structure described, in the local basis {eg = (1,0), ex = (0, -§z),e2 = (0, ■§-)}, as 



follows 

[(1, 0), (0, ^)] Rx tr 2 = (0, -F(^)), P. 0), (0, ^)1mxt R2 = (0, -^(^)). 

a a a a 

Prxtr 2 (1, 0) = 0, Prxtk 2 (0, g^) = Prxtr 2 (0, t^) = 

Therefore, = -F}, G 2 01 = -F 2 , & m = -J* C 2 2 = -F 2 , = 1 and p 2 = 1. 
Note that the homomorphism F generates a drag force F of the type 

F{x,y,A,y) = ^(x,y)— + F 2 (x,y) — ,F 2 1 (x,y)— + F 2 (x,y) — ^ 

= (Fl(x, y)(A -y)+ F 2 (x, y)(y + x), F$(x, y)(A - y) + F 2 (x, y)(y + x)) . 
On the dual bundle R x TR 2 we have a hamiltonian function: 

H(x,y,p x ,p y ) = -pi + -p 2 v +yp x - xp y + U{x,y) (4.2) 
and the corresponding Hamilton's equations are now: 

A = Px + y, 

y = Py~ X, 

OU OH 
Px = Vv--Q-~ F i ( x > y)Px ~ Ff (x, y)p y = -— F}(x, y)p x ~ F 2 (x, y)p v , 

OU OH 
P y = -p x ~ -q^ - F%(x, y)p x - F 2 (x, y)p v = -— - F 2 X (x, y)p x - F 2 (x, y)p y . 

Consider a section a € T(T*E 2 ) where a = a>i dx + a 2 dy. Then, 

dH d dH d 

c « = { Wx 0a) a^ + { dp- oa) dy 

Thus, Hamilton- Jacobi equation is equivalent to 
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For a hamiltonian function given by (4.2), the last two equations can be written as: 

dU , dai . . dcti . 

-^-a 2 + (a 1 +y)— + (a 2 -x)— + a i F 1 = 

dU . . da 2 , \ da 2 

— +a 1 + (a 1 +y)— + (a 2 -x)—+a l F 2 = 0. 

An interesting case is when the drag force is 

F(x,y,x,y) = (k(x,y)(x - y), k(x,y)(y + x)) , 

with fc G C*°°(]R 2 ). In this case, the homomorphism is F(X) = k(x,y)X with X G T( Xi? ,)M 2 . Thus, 
the equations of motion are 

x = p x + y, 

y = Py-x, 

- -*-f -«,,,)». 

In this particular case the linear almost Poisson tensor on K x T*R 2 is given by 

_ d d d d , , . d d , . ^ 9 9 

n Rx T*R2 = A q h A 3 1" K(sC) yjPcc ^— A 1- k(x, y)p y — A — 

dx dp x dy dp y dp Q dp x dp dp y 

where (po,x,y,p x ,p y ) are standard coordinates on K x T*R 2 . 

Now, if the function k is constant and we choose a section a — dS where S : K 2 — > K is an arbitrary 
function, the Hamilton- Jacobi equation is 

d(H o a) + kdS = 0, 

which is equivalent to the suggestive equation: 

kS + H o dS = constant 

or, in other words, 

kS(x, y) + H(x, y, — — , — — ) = constant. 
ox ay 

In particular, for the hamiltonian function given by (4.2 1, we obtain the following partial differential 
equation: 

, „. . 1 fdS\ 2 lfdS\ 2 dS dS TT , . 

kS{x ' y) + -2{d^) + 2{dy) +y dx~ -*9y+U(x,y)= constant. 

Note that, this equation is the Hamilton-Jacobi equation as stated in Corollary |3.3| for the cocycle 
(kS,dS) G C°°(R 2 ) x fi 1 ^ 2 ) ~ r((M x TR 2 )*) when we consider the skew-symmetric algebroid 
t : K x TM 2 -> M 2 

Example 4.5. Hamilton-Jacobi equation for a particle on a vertical cylinder in a uniform 
gravitational field with friction. As another example we consider a particle of mass m constrained 
to move on a cylinder of radius r in a uniform gravitational field of strength g and assume also the 
existence of a frictional force acting on the system. 

The Hamiltonian H : T*(R x S 1 ) -> R is: 

H(x,9,p x ,p e ) = + J + m^a; . 

Am Zmr z 

The frictional force is modeled in our setting by the homomorphism F : T(R xS 1 ) -> T(R xS 1 ) given 

by 

_, d , a 



A UNIFIED FRAMEWORK FOR MECHANICS 

with (Ki,K 2 ) e R 2 . 

The corresponding equations of motion are 
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mx 

,„2/i 



Px 

Pe 

Px = -mg - Kip x 
Pe = -K 2 pe ■ 

In this case, we may consider the skew-symmetric algebroid t:Rx T(R x S 1 ) — >■ R x S 1 associated 
with the above homomorphism F : T(M xS 1 )-) T(R x S 1 ) defined as in (12. 121) . For this skew- 
symmetric algebroid we have that = (1,0) G T(R x T*(R x S 1 )) = C°°(R x V) x fi^R x S 1 ) is a 
1-cocycle and V = $- 1 (0) = T(R x S 1 ). 

Let us consider a 1-form a £ T(T*(R x 5 1 )). If locally a is given by 

a = a x dx + agd9 

with a x , € C°°(R x S 1 ), then the local expression of the vector field on R x S 1 is 

Moreover, the Hamilton- Jacobi equation for a <E Sl^R x S 11 ) is 

i^da + d(H o a) + F* a = 0, 



a x d ag d 
m dx m 89 



where d is the standard differential and F*a is the pullback of a by F (see Example 2.5 1. In local 
coordinates this equation becomes 



a x da x 



ag da x 



+ mg + K\a x 



m dx mr 2 89 

a x dag ag dag 

7T~ H 2^0 + K2010 

m ox mr z 89 



In the particular case when a — dS with S a function given by S(x, 9) — (x) + (9), we have 
that the corresponding Hamilton- Jacobi equation is 



dx 



mg 



K. 



d9 



1 dSW d 2 SW 
m dx dx 2 

i ds (2 U 2 s^ 

mr 2 d9 d9 2 



Solving the equation we obtain that 

S {2 \9) = OorS (2) (0) 



K 2 mr 2 
2 ' 



gm 



—gm — gmW 

S^(x) = 

S m (x) = ±y / 2 v / -gm 2 x + C 2 if K i = 



if K x j= or 



where W is the Lambert W-function (the inverse function of /(«) = ve v ) and Ci and C 2 are constants. 

In Figure 1, we compare the trajectory of the particle for the free problem and the trajectory for 
the problem with friction. 
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Figure 1. Comparison of a free trajectory (without friction), on the left, and a 
trajectory with friction, on the right 



4.1.2. Unconstrained mechanical systems on Lie algebroids with a 1-cocycle. Let ([•,•], p) be a Lie 
algebroid structure (or more generally a skew-symmetric algebroid) on a vector bundle te : E — > Q 
and 4> £ T(E*) be a 1-cocycle such that <j>(q) ^ 0, for all q G Q. Denote by A (respectively, V) the 
afSne (respectively, vector) subbundlc of E given by A = </> _1 (l) (respectively, V = </> _1 (0)). 

In addition, we endow the vector bundle with a bundle metric S : E Xq E — >• R on £\ Denote by 
bg : — > E* the isomorphism of vector bundles induced by S- Consider the section X of E defined as 
follows 

X = bg 1 O (j). 

We will suppose, without loss of generality, that S(X, X) = 1. Thus, <p(X) — 1 and X is a section 
of the affine bundle tj\ : A — > Q. On the other hand, Q(X(q),v) = 0, for all v G V q , therefore 
E q =< X(q) > ®V q , for all q£Q. 

Now, let us consider the hamiltonian section h : V* -> _E* of the AV-bundle : -> V* 
characterized by 

h (Vq)(Vq) = Vq(Vq), h(jl q )(X q ) = ~H(r] q ), 

where H : V* — > M is the function 



V* -> E* of the AV-bundle fj, 
for all q £ Q, rj q G y g * and w 9 G V g , 



with SJ, : r x 7 



;S^(%,%)+V(q), 
the bundle metric induced by Q on V* and 



J^Ra real C°°-function on 
Q. Then, the function Fh ■ E* — > K associated with the section /i is just Fh — H o + X. 

Let (g l ) be a system of local coordinates for Q and consider an orthonormal local basis {eo, e a } of 
r(E) with e = X. Denote by (q l ,po,p a ) the local coordinates on E* with respect to the dual basis 
of {e ,e a }. 

The local expression of the hamiltonian section h G r(/i) is the following 
Kq\Pa) = {q\-H = - l -( Pa f -V(q),p a ). 

The integral curves of the hamiltonian vector field Rh G X(V*) are the solutions of the Hamilton 
equations 

= P0+ PaPa 



dt 

dpb = _ tdV 

dt Pb dq l 

For this mechanical system the dissipative term has the local expression 



{Ho[i,F h } = p l - 



GobPcPb- 
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Let a be a section of V * . Then, the section of E is given by 

Ch = iv o Ch + X, 

where Ch es the section of V defined by = P V (H) o a, for all /3 G r(V*). 



Thus, using Theorem 3.1 we deduce the following corollary 

Corollary 4.6. Let a be a section ofV*. Then, the following statements are equivalent: 

(i) If c : I — > Q is an integral curve of R% — Try* o R h o a G then a o c : / -> F' is an 

integral curve of Rh & X(V*). 
(ii) a € r(V*) satisfies the Hamilton- Jacobi equation: 

i^d v a + /jo ixd E (h o a) = 0. 



Remark 4.7. If /3 is a section of £*, the Hamilton- Jacobi equation for /x o f) is (see Corollary 3.2) 

i c »opd v (n qj3)+/jo ixrf B (/3) + o/io/3 + Xo/3)=0. 



If /3 is 1-cocycle on from Corollary |3.3| then the Hamilton- Jacobi equation is 

d v (H o no ft + Xo (3) = 0. (4.3) 



Therefore, Corollary 4.6 is a generalization of the main result of [28] (see Theorem 3 in 28\). In such a 
paper the authors obtain a Hamilton- Jacobi equation for mechanical systems on Lie affgebroids with 
this extra hypothesis on j3. 

If, additionally, V is a transitive Lie algebroid (that is, pv(Yq) — T q Q, for all q € Q) and Q is 
connected, we have that the Eq. (4.3) may be rewritten as follows 

Ho^loI3 + Xo/3 — constant. 



Example 4.8. Time dependent classical systems. Let it : Q — > K be a fibration on a manifold 
Q and 77 the 1-form on Q given by 77 = Tr*(dt), where t is the standard coordinate on K. Consider the 
standard Lie algebroid on TQ. Then 77 is a 1-cocycle for it and the afline bundle A = ?7 _1 (1) = {d£ 
TQ/n(v) = 1} Q may be identified with the 1-jet bundle of local sections of ir. Note that the 
associated vector bundle V = ?7 _1 (0) is just the vertical bundle of ir 

Vtt = {u£ TQ/r)(v) = 0}. 

Now, we take h : V*tt — > T*Q a hamiltonian section of fj, : T*Q — > y*^. If the local expression of 
h is 

h(t,q\ Pl ) = faq^-Hfaq^pi^pi) 
then the associated hamiltonian vector field Rh on V*tt is given by 

d dH d dH d 



1 dt ' dpi dq 1 dq l dpi 
Thus, the Hamilton equations are just the time dependent classical Hamilton equation for h 

dq l _ dH dp, _dH 
dt dpi ' dt dq l 

Now, consider a section a of the vector bundle V*tt —> Q. Then, = dFh is a vector field on Q 
defined by 

^ a )=i3 v (F fl )o| 1 oa, for f3en\E). 
The Hamilton- Jacobi equation is 

(i c *d(ho a))\ Vn = 0. 

In the case when a is a closed 1-form on Q the Hamilton- Jacobi equation may be rewritten as 

(d(F h oa)) lVlT = d v ^(F h o a )^0, 
i.e., F? — Fh ° a is constant on the fibers of tt. 
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Finally, we analyze the case when ir is trivial, that is, Q = R x P with P a connected manifold 
and 7r is the projection on the first factor. Then, Vtt =1x TP and the section h may be identified 
with a time dependent Hamiltonian function H : R x T* P — > R. If a = dW, with W:Q->Ka real 
function on Q, then 

dW 

(F h oa)(t,q) = — t + H(t,dW t (q)) 

with (t, q) elx P. Here Wt ■ P — > R is the real function defined by W t {p) = W(t,p). In this case the 
local expression of the Hamilton-Jacobi equation is 

dW dW 

— + H(t, q\ ~~) = constant on P, 
i.e., the time dependent classical Hamilton-Jacobi equation pQ. 

4.2. Nonholonomic mechanical systems with affine constraints. Let ([■,■], p) be a Lie alge- 
broid structure on a vector bundle te ■ E — > Q. 

A mechanical system subjected to affine nonholonomic constraints on E consists of 

(i) a vector subbundle txj : U — > Q of E, 

(ii) a bundle metric Q : E Xq E — > R on £, 

(iii) a function V : Q — > R on Q 

(iv) and a section X e T(E) such that 7(X a ) = 0, where ? : E = £7 © J/- 1 -> f7 is the orthogonal 
projector defined by the metric S- 

Then, one may consider an affine bundle tu ■ U — > Q, 

q e Q — ► U q = {X (q) + u q /u q e U q } 

whose associated vector bundle is just U, describes the affine nonholonomic constraints. Denote U + 
the affine dual bundle associated to II whose fiber at q € Q consists in the affine functions over U q . 
Moreover, U + has a distinguished section <fi '■ Q — > which is induced by the constant function 
4> q — 1 on lig. 

On the other hand, if we denote by It = (U + )* the bidual bundle of U, then U is a vector subbundle 
of K x E -> Q with fiber at q £ Q 

U q = {(X, XX (q) + u q ) / X E R and u q E U q }. 

Thus, a section of U may be identified with a pair (f,fXo + a), with a € T(U) and / e C°°(Q). 
Under these identifications, the distinguished section <p is given by 

<f>{f, fX a +a) = f. 

Moreover, in a natural way, the projector 7 : E = U © U 1 - — > U defined by the metric S induces a 
new morphism 7 : R x E — > U of vector bundles given by 

7(X,e q ) = (X,XX + 7(e q )) 1 

for all A G R and e q e E q . 

In what follows, we will introduce a skew-symmetric algebroid structure on li such that is a 
1-cocycle. In order to do this, we consider the Lie algebroid structure ([•, -]r x _e, Prxe) on IR x E 
induced by the Lie algebroid structure on E and the homomorphism F = on E (see Example 2.5). 
Now, we consider the bracket [•, on the space of sections of U and the vector bundle morphism 
/Of; : U — > TQ given by 

KfiJiXo + ax), (/ a , f 2 X + a 2 )] Q = f ± X + a,), (/ 2 , / 2 X + a 2 )j MxE ) 

fX +a) = PuMf, fXo + <r) = p{f*o + a) 
for a,ai,a 2 E T(U) and f,fi,f 2 ,f € C°°(Q). Then, using that T:lxi?->[/isa projector, we 
deduce that the pair ([•, ^^Pu) 1S a skew-symmetric algebroid structure on IX. With respect to this 
structure, one may prove that 

d^cb = 0. 
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Note that the corresponding skew-symmetric algebroid structure on _1 (O) = U is just 
[o-i,o- 2 ]£/ = y([o-i,o- 2 ]), pu{e)=p(o-), with a t ,a G r(t/)- 
Moreover, CP : _E — > [/ and 3 : R x E 1 — >• It are skew-symmetric algebroid morphisms. 
On the other hand, one may consider the hamiltonian section h : U* —> \L + defined by 

h (v q )(^-AX (q) +u q ) = r)g(u q ) - XH(r) q ), Vr? g G U*, and (\,AX (q) + u q ) e U q , 
where H :[/*—> K is the real function 

#fa 9 ) = ^Sir.(»fo,»fo)+V(g). 

Here, Si/* : U* x U* — > R is the fiber metric induced by S on U*. In this case, we have that 
-F/j = (1, Xo) +H ofi, where (1, Xq) is the linear function on tt + induced by the section (1, Xo) £ T(IV). 

Let (q l ) be a system of local coordinates for Q and consider an orthonormal local basis {e a ,eA} 
of T(E) adapted to the decomposition E — U © U^. Then, {(1, Xq), (0, e )} is a local basis of 
sections of It. Denote by (q z ,Po,p a ) (respectively, (q l ,p a )) the corresponding local coordinates on lt + 
(respectively, U*) with respect to the dual basis of {(1, Xo), (0,e a )} (respectively {e a }). 

The local expression of the hamiltonian section h G r(/i) is the following 
h{q\p a ) = (q\-H = - l -( Pa f -Y(q), Pa ). 

The integral curves of the hamiltonian vector field Rh € X(U*) are the solutions of the Hamilton 
equations 

~dt = Po+PaPa 

where T([Xo,e;,]) = CQ b e c , p(X ) = Po^J" an d ^afc anc ^ Pa are local structure functions of E. A 
Lagrangian version of these equations was considered in [TS] . 

Now, let a be a section of U* . In this case, we have that the section of It defined as in (3.5 1 is just 

C,? = (i,G + x ), 

where is the section of U given by 

v((h) = v v (h)°u, V V eT(U*). 

Therefore, from Theorem |3.1[ we deduce that 

Corollary 4.9. For a G T(U*), the following statements are equivalent: 

(i) If c : I — » Q is an integral curve of R^ = Ttjj* o R h o a G £(Q), </ien a o c : I — > U* is a 

solution of the Hamilton equations, 
(ii) a G T(U*) satisfies the Hamilton- Jacobi equation: 



i(f,d u a + fio i( 1Xo )d u (h o a) = 0. 



Remark 4.10. The section ujf t = p, o in Xo \d(h o a) on U* can be written as 
u$(X) = p(X )(a(X)) + p(X)(H o a) - a(7(lX , X})). 

o 

From Corollary |3.3| we conclude that 

Corollary 4.11. Suppose that Q is a connected manifold and that the finitely generated distribution 
V defined by V q := pu(U q ) for all q G Q, is a completely nonholonomic distribution. If fi is a section 

of IL + such that d u fj = 0, then the following statements are equivalent: 
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(i) If c : I — > Q is an integral curve of Rf°^ = Ttjj* ° Rh ° /J ° P & then /io/ioc:/- > V* 

is a solution of the Hamilton equations, 
(ii) /3 £ r(U + ) satisfies the Hamilton- Jacobi equation: 

H o /! o j3 + /9(1, Xq) = constant. 

Example 4.12. An homogeneous rolling ball without sliding on a rotating table with time- 
dependent angular velocity. We consider a homogeneous ball with radius r > 0, mass m and 
inertia mk 2 about any axis. Suppose that the ball rolls without sliding on a horizontal table which 
rotates with a time-dependent angular velocity tt(t) about vertical axis through of one of its points. 
Apart from the gravitational force, no other external forces are assumed. 

Choose a cartesian reference frame with origin at the center of rotation of the table and z— axis along 
the rotation axis. If (t, q 1 , q 2 , q 1 , q 2 , u>i, U!2, uj 3 ) are the corresponding coordinates over R x TR 2 x R 3 , 
then (<? 1 ,<; 2 ) denotes the position of the point of contact of the sphere with the table and o>i,W2 and 
W3 are the components of the angular velocity of the sphere. 

Note that since the ball is rolling without sliding, then the system is subjected to the affine con- 
straints 

q 1 - ruj 2 = -n(t)q 2 



+ rw 1 = iltyq 1 . 



Let (t, q 1 , q 2 ,Pi,P2, ^1,^2, ^3) be the corresponding coordinates on (MxTIR 2 xIR 3 )*. The hamiltonian 
section h : (K x TM 2 xR 3 )MIx(lx TM 2 x I 3 )* of the system is given by 

h(t,q 1 ,q 2 ,Pi,P2,Tr 1 ,ir 2 ,Tr s ) = (-H(t,q 1 ,q 2 ,p 1 ,p 2 ,Tr 1 ,TT2,7T 3 ),t,q 1 ,q 2 ,p 1 ,p 2 ,Tr 1 ,TT2,TT 3 ) 

where H : (E x TR 2 x R 3 )* ->• R is the real function 

h = \{-(pI+p 2 2) + \M + A + 

l m mk z 
Moreover, the constraints may be rewritten as follows 

1 r 

ipl = il{t)q 2 + —pi — 7T 2 = 

m mk 

1 r 

i>2 = -flitjq 1 + —P2 + —p;TTl = 0. 

to mk z 

Then the Hamilton equations of this non-holonomic system are 

1 1 

q = —Pi 

777 

2 1 

q = —P2 



m 

mk 2 dVL{t) 2 n , ,l>2. 



mk 2 dCljt) ! -.ft. 

P2 = WTr- 2i ^r q +m m ] 

rmk 2 dVt{t) 1 nM Pi. 
-1 - ^2(^9 +«(*)-) 

rmk 2 dVt{t) 2 nM P2, 

P3 = 

and the constraints =tp2 = Q (for more details, see [18]; see also [38]). 

Our goal is to encode all this information in a mechanical system subjected to affine nonholonomic 
constraints on a Lie algebroid. Consider the vector bundle te ■ E —> Q, where E :- 
and te '■ E — > Q is defined by 

TE(t,q 1 ,q 2 ,i,q L ) q 2 ,ui,u>2,va,) = (t,q X ,q 2 )- 
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We choose the following global basis of T(E) 

e = (l-m^ + m^0) ^ = (^,0), e 2 = (^0), 

e 3 = (0,(1, 0,0)), e 4 = (0,(0, 1,0)), e 5 = (0,(0,0,1)), 

On E we define the Lie algebroid structure by 

[e ,ei] = -ft(i)e 2 , [e , e 2 ] = to(t)ei, [e 3 , e 4 ] £ = e 5 , [e 4 , e 5 j E = e 3 , [e 5 , e 3 ] £ = e 4 , 

d d d d d 

p E (e ) = -Q t -n(t)q 2 — + n(t)q 1 — , p B ( ei ) = — , p £ (e 2 ) = — . 

The rest of the local structure functions are zero. 
The constraints induce a vector subbundlc of E 

U := span{e 3 — re 2 , e 4 + re\ , e 5 }. 

Consider the bundle metric on E 

S = e 2 + M( ei ) 2 + (e 2 ) 2 ) + mfc 2 ((e 3 ) 2 + (e 4 ) 2 + (e 5 ) 2 ). 

In order to do the decomposition E = U ® t/- 1 , we take the following orthonormal basis of T(E) 
with respect to 9 

e"o = e , ei = (re 3 + fc 2 e 2 ), e 2 = (re 4 - fc 2 ei), 

1 1 / x - 1 

e 3 = — — (e 3 — re 2 ), e 4 = — — (e 4 + rei), es = - — ^es. 

Vm(fc 2 +r 2 ) V v / m ( fc2 +^ 2 ) fc ^ 

Then, {e 3 ,e 4 ,e 5 } (respectively, {e ,ei,e 2 }) is a orthonormal basis of U (respectively, U^). 

Moreover, for this mechanical system, the distinguished section X of E is X = eo- Note that 
3>(X ) - 0. 

Denote by (tjq 1 ,q 2 ,po 1 pi 1 p2,Tti, ^2,^3) the coordinates on E* with respect to the dual basis 
{e ,^ 1 ,^ 2 , e 3 ,e 4 ,e 5 } of {e , ei, e 2 , e 3 , e 4 , e 5 }. With respect to these coordinates the function : 
U* -)• R is defined by 

if(7fl,7f 2 ,7T 3 ) = -(7f 2 +7T2 +^ 3 ) 

and the structure functions of the skew-symmetric algebroid on It — > R 3 with respect to the basis 
{(1,X Q ), (0,ei)}i = 3,4,5 are the following 

k 

L 34 -L 45 -L 53 - ^_ (r2 + p) , 



p8 = i, - -o(t)g 2 , p 2 = n(*y, pi = -p§ = 



yfm(r 2 + k 2 ) 

Let us consider the section a £ T(U*) to be a — d u g for the real function on R 3 

9 = gfaq 1 ^ 2 ) = <Pi(t)<l 1 + <P2(t)q 2 

where (fii,ip 2 e C°°(R). Then, 

y/m(k 2 + r 2 y dq 2 dq 1 1 ^m(k 2 +r 2 Y * V ' } 
and the section e r([/) is 

Ch = , „ 2 , 2 , (-y 2 (t)e 3 + ^i(t)e 4 ). 
\/m\k 2 + r z ) 

It is important to note that a G L(J7*) is not a 1-cocycle of the skew-symmetric algebroid tu* : 
U* Q. In fact, 



24 



P. BALSEIRO, J. C. MARRERO, D. MARTIN DE DIEGO, AND E. PADRON 



However, 



ic* d u a = 0. 



Thus, the Hamilton- Jacobi equation becomes 



n{t)r 2 

k 2 + r' 



■Mi)- 



(4.4) 



Now, in order to apply Corollary 4.9 we have to find an integral curve c(s) = (t(s),q 1 (s),q 2 (s)), 
for s £ R of the vector field i?£ € X(Q) given by 



d 



d 



dq 1 \m(k 2 + r 2 ) 



dt \m(k 2 + r 2 ) 
Then, the curve c has to verify t(s) = s + cq, and taking Co = we get 

_2 



if 2 + ttitfq 1 



d 

dq^' 



fit) 



m(k 2 + r 2 ) 



m(k 2 + r 2 ) 



We conclude that aoc(t) = (t, q l {t) 1 q 2 (t); 



Vl {t)~n{t)q 2 {t) 
ip 2 (t) + n(t)q 1 (t). 



y/ m(k 2 + r 2 ) \J m(k 2 + r 2 ) 



(4.5) 

(pi(t),0) is an integral 



curve of Rh, where ifi(t) and q l (t) are real functions that satisfy (4.4) and (4.5). 

As a particular case, we can take the angular velocity of the table to be ft(t) = Qq = cte > and 
we get that the curve a o c(t) = (t, g 1 (t), q 2 (t); at3(c(t)), Q!4(c(t)), 0) is given by 



a 3 (c(t)) 



a 4 (c(t)) = 



a/ m(k 2 + r 2 ) 
r 

yjm{k 2 + r 2 ) 



C\ sin 



C\ cos 



r 2 n Q t 

k 2 + r 2 

r 2 n t 



C2 cos 



— Ci sin 



r 2 n t 

k 2 + r 2 

r 2 Q t 
k 2 + r 2 



and g 1 (t),g 2 (i) solutions of the system (4.5). The trajectories of the ball on the rotating table (tra- 
jectories in (qi(t),q2(t))) arc ellipses centered in the origin of the table, which depend on the initial 
conditions of the problem. 

If fl(t) = Clot then the curve a o c(t) = (t, q 2 (t); as(c(t)), 0:4 (c(i)), 0) is given by 

r 2 n t 2 \ , „ ( r 2 Q t 2 



a 3 (c(t)) = 



a 4 (c(t)) 



v /m(k 2 + r 2 ) 
r 

\J m(k 2 + r 2 ) 



C\ sin 



Ci cos 



2(k 2 + r 2 ) 

r 2 fl t 2 
2{k 2 + r 2 ) 



C2 cos 



\2(k 2 



C 2 sin( 



2(k 2 + r 2 ) 



where C\,C2 are real constants. In this case, the solutions (qi(t), 52(i)) of (4.5) give trajectories on 
the table as in Figure 2. 



FIGURE 2. The trajectory of the ball on the plane with velocity Q(t) = t for < t < 5 
and < t < 30. In the last two figures the velocity is changed Q(t) — lOt for < t < 20 
and < t < 30 
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4.3. A example of a nonholonomic mechanical system with linear external forces: the 
vertical rolling disk with external forces. We will use the classical example of the vertical 
rolling disk to show how external forces can be encoded in the geometric structure of the constraint 
submanifold. Then, we are going to find the Hamilton-Jacobi equation and we obtain some particular 
solutions. 

Consider a vertical disk that is allowed to roll on the xy-plane and to rotate about its vertical 
axis. Let x, y denote the position of contact of the disk with the xy-plane, 9 will denote the rotation 
angle of a chosen point P of the disk with respect to the vertical axis and finally 4> will represent the 
orientation angle of the disk as in Figure 3. 



Therefore, the configuration space for the rolling disk is Q — M 2 x S 1 x S 1 with coordinates 
(x,y,9,<p). On the tangent bundle TQ — > Q we consider the Lie algebroid structure ([•, -]tq, idro) 
where [•, -}tq is the usual Lie bracket of vector fields and the anchor map, in this case, is idrQ- 

The Lagrangian for this system is: 



where m is the mass of the disk, / its moment of inertia about the axis perpendicular to the plane 
containing the disk and J is the moment of inertia about an axis in the plane of the disk. This 
Lagrangian induces a fiber metric 



The nonholonomic constraints of rolling without slipping are 

( X = (i?COS <f)6 

{ y = (R sin <f)6 

and they define the constraint subbundle td : D — >• Q of TQ. 

In terms of the fiber metric 9, we find an adapted basis for the nonholonomic problem. More 
precisely, we look for an orthonormal basis of vector fields {X 1} X 2 , X s , X 4 } of TQ such that D — 
span{X i,X2} and D^- = span{X 3 ,X4}. This basis is given by 




v 



Figure 3. The vertical rolling disk 




3 = m(dx ®dx + dy® dy) + Id9 ® dQ + Jdcj) <g> d(j>. 




X 2 



J_d_ 




26 P. BALSEIRO, J. C. MARRERO, D. MARTIN DE DIEGO, AND E. PADRON 

We endow the fiber bundle t& : D — > Q with a skew-symmetric algebroid structure ([•, Jd, Ad) 
defined by (see Example 2.6 ) 

\X l ,X 2 \ D =P{[X u X 2 ] TQ ) and p D (X x ) = X 1; p D (X 2 ) = X 2 , 

where P : TQ — > D is the orthogonal projector (with respect to the decomposition TQ = D D- 1 ). 
Note that, po = Ptq ° *d with : D TQ the natural inclusion. Therefore, in terms of the basis 
{Xi, X 2 }, the (non zero) local structure functions of the skew-symmetric algebroid on D are given by 

Rcoscj) e 1 

\'mR z + 7 \/rriR z + 1 

y i?sin0 ^ 1 

(™ = / p2 , r ' ^ D ^2 = -7f • ( 4 -6) 

Since [Xl, X 2 ]tq G span{X 3 } we have that Q{ 2 = Q 2 2 = 0- 

In coordinates (v 1 , v 2 , u 3 , v 4 ) induced by the orthonormal basis of sections {Xi, X 2 ,X%, X4} the 
Lagrangian is 

L(x, y, 9, 0; v\v\v\ v 4 ) = \ ((^) 2 + (v 2 ) 2 + (v 3 ) 2 + (« 4 ) 2 ) , 

and the equations determining the constraints are v 3 = v 4 — 0. Therefore, the restricted lagrangian 
L D : D -> R becomes L D (x, y, 0, 0; u\ v 2 ) = \ ((v 1 ) 2 + (v 2 ) 2 ) . 

Consider now, the dual vector bundle m* ■ D* — > Q with coordinates (a:, y, 0, 0;pi,P2) induced 
by the dual basis {X^X 2 } of {Xi,X2}. Then, the vector bundle td» : D* —> Q has a linear almost 
Poisson structure given by 

i?cos</> 1 



^mR 2 + r ^JmR 2 + f 

Rsincj) 1 

{y,Pl}D* = —j====, \<P,P2)D- = —j= 

and the other fundamental brackets are zero. 

In these coordinates, the Hamiltonian function H : D* M can be written as 

1 



H(x,y,8,(t>;pi,p 2 ) = - ((pi) 2 + {p 



2 , S)2) 



It is very interesting the study of the rolling disk with external forces [35]. The system has two 
natural inputs, a torque that makes the disk spin and another one that makes the disk roll. First we 
are going to study the most general situation and then we will analyze particular cases. Suppose that a 
linear force is acting on the disk, then the pull back of this force in D* is given by F(q, v) — (F 1 (q)v 1 + 
F^(q)v 2 )X 1 (q) + +F 2 (q)v 2 )X 2 (q), where (g,t>) = (x, y, 9, <j>; v 1 , v 2 ) and F? £ C°°(Q). 

Since the chosen basis is orthonormal, we have that the homomorphism F : D — » D induced by the 
force F is 

FiX,) = F^X 1 +F 2 X 2 
F(X 2 ) = F 2 1 X 1 +F 2 X 2 

and thus the skew-symmetric algebroid onlxD has (non z ero) local structure functions given by 
Gj x = -Fl, C 2 ! = -F 2 , & 02 = -F 2 \ C 2 2 = -F 2 and equation @. 

Therefore, the corresponding Hamilton equations modified by the action of an external force are 

R cos <j> R sin <p 

Pi, y= r-f — =^PU 



y/I + mR 2 ' VI + mR 2 

1 X 1 



Vl + mR 2 y/1 
Pi = -Flpi - F 2 p 2 p 2 = -F 2 pi - F 2 p 2 . 



In order to write the Hamilton- Jacobi equations, let us consider a section a € T(D*). 



A UNIFIED FRAMEWORK FOR MECHANICS 



27 



Then, such equations are 

a\.X\{ai) + a 2 .Xi (a 2 ) + a 2 X 2 (a 1 ) — a 2 .X\ (a 2 ) + a.\.F\ + a 2 .F 2 = 



ct\.X 2 (a\) + a 2 X 2 (a 2 ) — a 1 X 2 (a 1 ) + a 1 X i (a 2 ) + a\.F 2 + ol 2 .F,^ 







(4.7) 



where a = u\X x + a 2 X 2 and cei, a 2 £ C°°(Q). 

Particular case: A torque that makes the disk spin. 

Let us consider the external force F = X(4>)4>d(j), with A £ C°°(R). Writing this force in terms of 
the dual basis {X X ,X 2 } we obtain 

F( q ,v) = X fv>X>, 
where (q, v) — (x, y, 0, <f>, v 1 , v 2 ). Therefore, the homomorphism F : D — > D is 

F(X 1 )=0 and F(X 2 ) = ^X 2 

and the ske w-sy mmetric algebroid onKxD has (non zero) local structure functions given by C§ 2 = 
-^f - and (46 1. 

Consider a section a £ T(D*) such that a = kX 1 + a 2 (4>)X 2 , with k = constant. 
Thus, Hamilton- Jacobi equation (4.7) is simply (note that, in this case, d E a = 0), 



Therefore, from (4.8), we deduce that 

a 2 {4>) = 

where k is an arbitrary constant. 



K<t>) 
VI' 



X(s)ds + k 



(4.8) 



By Eq. (3.1), we have 

R k cos 4> d Rk sin 4> d 

Rh = 



d 1 



h VI + mR? dx VI + mR 2 dy Vl + mR 2 d6 J 
We conclude, by Corollary |4.1| that 



X(s)ds — k 



d_ 
d4>' 



aoc{t) = {x(t),y(t),9{t),ct>{t);k,- 



'J Jo 



X(s)ds — k) 



is an integral curve of Rh £ X(D*), if e(t) = (x(t),y(t),6(t),<p(t)) is an integral curve of R^. 

As a particular case, we fix X(<p) — Kcostfi with K = cte ^ 0. Hence by equation (4.8) we 
have that a 2 (cf>) = —-^=sixi(j) + K but, just for simplicity, we will choose k = 0. If c : I — > Q, 

K 

c{t) = (x(t),y(t). z(t),6(t),(j)(t)), is an integral curve of R^ then 4>{t) = — — sine/). That is, 



<j>(t) = 2arctan fe~* i+< M 



with (j)Q an arbitrary constant. Therefore, the solution of the system, modified by an external force 
F = (K cos (f>)(j> d<f> that makes the disk spin, is 

ao C (f) = {x(t),y(t),9(t),<f>{t);k,---=sm<f)(t)), 
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where x(t),y(t),Q(t),<f>(t) are curves given by 

*> " 7TTSJp(' + i ta ( 1 + «- ,Wfc )) + - 

= 2 arctan(e _ 

where xo, yo, 9q, <f>o are arbitrary constants. 

We also have the dissipative term for this case given by 

{if o^Fj = ^(Par- 
Remark 4.13. The function / <= C°°(Q), given by 

><«=i?(f A(,)<! *) = -^ sinV 

verifies that F*a = d D f. Thus, we obtain that the Hamilton Jacobi equation can be written as 

K 2 

H o a — — — sin 2 <f> = constant, 
on Q, since D is a completely nonholonomic distribution. 



5. Conclusions and future work 



A Hamilton-Jacobi equation for a great variety of mechanical systems is derived. The type of 
systems considered includes mechanical systems with dissipative forces, nonholonomic system sub- 
jected to linear or affine constraints or, even, explicitly time-dependent mechanical systems. With 
this general purpose in mind, we find that the geometric structure of skew-symmetric algebroid has 
the appropriate inclusive nature, adequate to model all these different types of mechanical systems. 
Adopting this point of view we prove a general version of the Hamilton-Jacobi equation for skew- 
symmetric algebroids with a distinguished cocycle, specializing the results for the different mechanical 
systems under study. Several examples prove the utility and novelty of our results. 

Of course, a lot of work must be done in future research. For instance, in our paper a crucial 
assumption is made: all the constraints are linear or affine, even the dissipative forces considered 
are of a very special type (in such a way that they induce a linear bivector on the dual bundle). It 
would be interesting to discuss the more general case in a non-linear setting, discovering the underlying 
geometric structures and deriving, if possible, a Hamilton-Jacobi equation. Moreover, in future papers, 
we will study more explicit examples of applications of our theoretical setting, analyzing when the 
separation of variables technique works and relating it with topics like integrability. Also, our setting 
is ready for the introduction of control forces and therefore for the study of controlled mechanical 
systems and, as a consequence, to address problems like kinematic reduction, kinematic controllability, 
Hamilton-Jacobi-Bcllman equation in optimal control, etc. 
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